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ABSTRACT 

We reanalyze the catalogs of molecular clouds in the Local Group to determine the parameters 
of their mass distributions in a uniform manner. The analysis uses the error-in-variables method of 
parameter estimation which accounts not only for the variance of the sample when drawn from a parent 
distribution but also for errors in the mass measurements. Testing the method shows that it recovers 
the underlying properties of cumulative mass distribution without bias while accurately reflecting 
uncertainties in the parameters. Clouds in the inner disk of the Milky Way follow a truncated power- 
law distribution with index 7 — — 1.5±0.1 and maximum mass of 10®'^ Mq. The distributions of cloud 
mass for the outer Milky Way and M33 show significantly steeper indices (70 mw — —2.1 ± 0.2 and 
7M33 = — 2.9±0.4) with no evidence of a cutoff. The mass distribution of clouds in the Large Magellanic 
Cloud has a marginally steeper distribution than the inner disk of the Milky Way (7 = — 1.7± 0.2) 
and also shows evidence of a truncation with a maximum mass of 10^'^ Mq. The mass distributions of 
molecular clouds vary dramatically across the Local Group, even after accounting for the systematic 
errors that arise in comparing heterogeneous data and catalogs. These differences should be accounted 
for in studies that aim to reproduce the molecular cloud mass distributions or in studies that use the 
mass spectrum as a parameter in a model. 
Subject headings: ISMiclouds — methods:data analysis 



1. INTRODUCTION 

The mass distribution of molecular clouds is one of the 
primary characteristics of the their population. In the 
inner disk of the Milky Way, the mass distribution fol- 
lows a power law with dN oc M'^dM, 7 ~ —1.5. More 
recent surveys of molecular clouds throughout the Local 
Group find that the mass spectrum also follows a power- 
law, but t he indices a rc steeper than in inner Milk y 
Way (e.g. lEngarglola et al.ll2003t iMizuno et al.ll200H) . 
Indeed, the mass spectrum may be the only feature 
of the molecular cloud population that varies between 
systems, since other cloud properties (e.g. cloud radius 
and line width) obey th e relationship s established in the 
Milky Way fe lson & ScoviUe" 1993 IR,osolowskv''eFaLl 
120031 IMizuno et al.ii2001i) . Careful attention to accu- 
rately determining the parameters of the mass spec- 
trum is critical in using the mass spectrum to quan- 
tify differences between cloud populations. In addition, 
the empirically derived mass distribution is an impor- 
tant parameter for theoretical and modeling work. Sev- 
eral studies aim t o reproduce t he mass distribution of 
molecular clouds fKwan| | 197£ jElmegreen & Falgarond 
^996; Vazaucz-Semade ni et al.l ll997: Stutzki et al. 199^^ 
IWada et al.l |2(}00) or use the mass spect r um as in- 
puts to models LMcKee fc Wilhamsl Il997t iTanI 120001 
iKrumholz fc McKed \200^\ . Most of these studies fo- 
cus on the canonical value of 7 « —1.5 adopted from 
the inner Milky Way, neglecting any variation in the dis- 
tribution. Judging from the scope of these other stud- 
ies, measuring the mass distribution of molecular clouds 
is essential for understanding both cloud formation and 
the importance of star-forming clouds in regulating large 
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scale star formation. 

Since the parameters of the cloud mass distribution 
are widely used in the study of the star-forming inter- 
stellar medium, this paper outlines some of the pitfalls 
associated with the standard methods of estimating the 
parameters of power-law distributions and suggests im- 
provements to minimize inaccuracy O. With these im- 
provements, we reanalyze data from existing catalogs of 
molecular clouds (33) and note interesting results {^^■ 
This work stresses the importance of accounting for the 
observational uncertainties and systematic effects that 
bedevil the study of molecular clouds. Accurately deriv- 
ing the index of a power-law distribution is also useful 
for studying populations of other objects. In particu- 
lar, the derived mass spectrum of clumps within molecu- 
lar clouds is subject to identical systematics as the mass 
distributions studied in this work. The methods devel- 
oped in this study as well as their attendant cautions are 
directly applicable to the study of clump mass distribu- 
tions and their relevance in the formation on individua l 
stars (e.g. Will iams et a l. 1994: Stutzki fc Giist en"l990'). 
Luminosity and mass distributions of stars and galax- 
ies are characterized by non-linear distributions and the 
techniques presented in this paper readily extend to the 
study of these objects. 



2. FITTING MASS SPECTRA 

The mass distribution of a population of molecular 
clouds is usually expressed in differential form, namely 
the number of clouds that would be found in a range of 
masses. In the limit of a small mass bin, this is expressed 
as 



dM ' 



(1) 
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This expression can be integrated to give the cumulative 
mass distribution 

N{M' > M) = / f{M)dM = g{M), (2) 

which gives the number of clouds with masses greater 
than a reference mass as a function of that reference 
mass. For molecular clouds, both forms of the mass spec- 
trum obey power-laws: f{M) (x and g{M) oc M'''+^ 
with 7 < —1 in all known cases. Some mass distribu- 
tions lack clouds above some maximum mass Mq. To 
account for this feature, we a dopt a truncated power-law 
distribution as suggested bv IWilliams fc McKed l|1997fl 
and alter their formalism to our notation. The full form 
of the cumulative distribution is 



N{M' > M) ^ No 
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where Mq the maximum mass in the distribution. Nq is 
the number of clouds more massive than 2i/(T+i)Mo, the 
point where the distribution shows a significant deviation 
from a power law. If A^o ~ 1, there is no such deviation. 
For this form of the cumulative mass distribution. 
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In most studies, only the index 7 is reported since A^o is 
assumed to be 1 and Mq is the maximum mass cloud in 
the sample. The index is the most important parameter 
since it describes how the integrated mass is distributed 
between the high and low mass members of the cloud 
population. For values of 7 > —2, the majority of the 
mass is contained in the high mass clouds and the reverse 
is true for 7 < —2. When 7 < — 2, the integrated mass di- 
verges as M — ^ 0, implying a break in the power-law be- 
havior of the mass spectrum at or below the completeness 
limit to ensure a finite integrated mass. Distributions 
with A^'o > 1 are also physically interesting since they 
have a characteristic feature in an otherwise featureless 
mass d istribution. In the Milky Way, Williams & McKee 
l)1997|) report evidence that Nq is significantly different 
from unity, implying a cutoff at high mass (3 x 10® Mq) 
in the Galaxy. The parameters of the mass distribution 
are important both as predictions of theories as well as 
inputs to models. It is critical to estimate these parame- 
ters with minimum bias from the mass measurements of 
a cloud population. 

2.1. Binned Mass Spectra 

Most studies of the mass spectrum of giant molecular 
clouds (GMCs) estimate the slope of the mass spectrum 
by fitting an approximation of the differential version. 
They generate this approximation by separating the mass 
measurements into logarithmically spaced bins. Then, 
the number in each bin {Nf,in) is divided by the width of 
the bin AM: dN/dM « Nun/AM. The uncertainties 
in these bins are then assumed to arise from counting 
errors, so ahi ri = \/Nhin.l AM. Studies using this tech- 
niquc include 'Solomon ct al." ('198/): 'Wil liams fc McKed 
^997); Engargiola ct al. (2003); Mizmi oeTalJ l|200lD : 
iHever et al.l (j2001f) and many others. 

There are two principal drawbacks to this technique: 
(1) it is sensitive to the selected values of bin size and 



bin spacing and (2) it neglects errors in the mass deter- 
mination of the clouds, which can be substantial. Figure 
^ shows the variation in the derived index of the mass 
spectrum for different choices of bin size and bin posi- 
tion. To generate t hese figures, we used the mass data 
from iSolomon et al.i (JDS?, SRBY) with the same com- 
pleteness limit of 7 X IO^Mq as is quoted in their paper. 
For a given set of bin parameters, we fit a power-law dif- 
ferential mass spectrum to the results to all data that are 
at least one full bin above the completeness limit. We fol- 
low the method of Williams fc Mc Kcc ( 1997) for the fit 
and the determination of errors in the mass distribution. 
The systematic error in the parameters is comparable the 
errors typically quoted in these studies. Such errors be- 
come negligible in the limit of large numbers of clouds. 
In the study of He ver et al.l ( 2001, HCS), there are over 
1300 clouds above the completeness limit as opposed to 
only 200 in the SRBY study. When the same experiment 
is performed on this much larger sample, the variation in 
the derived index reduces to ±0.05 and agrees with the 
— 1.8 quoted in the HCS paper. To use binned mass spec- 
tra in estimating the parameters of the mass distribution, 
the sample should have Ndouds > 500 to reduce errors 
to less than 0.1 in the index. 

In addition to large variations in derived bin param- 
eters, the binning method also neglects the principal 
source of uncertainty, namely the mass measurement it- 
self. The mass of a molecular cloud is notoriously dif- 
ficult to calculate. The principal methods for deriv- 
ing the mass are using the CO-to-H2 conversion fac- 
tor and using the virial theorem. The conversion fac- 
tor linearly scales the integrated CO surface brightness 
to a column density along a line of sight. With a dis- 
tance measurement, the column density of the cloud 
and the area on the sky are combined to calculate the 
cloud mass. The conversion factor is empirically tested 
to trac e H2 column density across a variety of environ- 
ments (jBloemen et al.| 1986ft t hough variation is reported 
among galaxies (jArimoto et al. 1996). Within a single 
galaxy, however, the conversion factor has been found 
to be constant despite c hanges in the galactic environ- 
ment (|Rosolowskv et al.ll2003() . The virial mass mea- 
surement assumes that clouds are virialized and uses the 
resolved cloud sizes {Re) and line widths (ay) to convert 
to a virial mass: Myx — 5Re(Jy/{aG) where a is the 
virial parameter, which depends on the mass distribu- 
tion within the cloud as well as the influence of magnetic 
fields and external pressure on the energy balance of the 
cloud. HCS present evidence that molecular clouds with 
Al < 10* Mq are not virialized and the virial mass mea- 
surement overestimates the masses of these clouds. 

Both of these methods for measuring cloud masses are 
subject to large (< 50%) errors. Absolute flux cali- 
bration of CO data is rarely accurate to better than 
10%. The variations of the conversion factor with 
physical conditions remain poorly understood, despite 
many attempts to q uantify them fWol fire et alJ Il993t 
IDickman et aljri986f) . Finally, the distances to most 
molecular clouds are difficult to measure. For Milky 
Way molecular clouds, most distances are determined 
kinematically with the distance degeneracy for the in- 
ner Galaxy being broken by angular scale, displacement 
above the plane, and association with other objects of 
known distance (SRBY). Distance measurements are also 
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Fig. 1. — Demonstration that fitting binned mass speetra is sensitive to clioice of bin size and offset. The left panel shows the derived 
index of the mass spectrum for the clouds in Solomon et al. ( 1987) as a function of the width of the logarithmic bins (in dex) used in fitting. 
The fits used the same completeness limit as the original study (7 X IO^Mq). The bin size ranged from a quarter octave 0.25 log ]^q(2), to 
a full octave, logj^Q(2). The right panel shows the variation in mass spectrum for half octave bins and different bin positions. The bins are 
shifted in log-space by the quoted bin offset. Both variations in parameters show significant variation in the derived index. 



important in measuring virial masses since the physical 
size of the cloud is determined by converting an angular 
scale to a physical length. In addition, small or distant 
clouds are often poorly resolved and great care must be 
taken to measure the radius of an intensity distribution 
that has been convolved with the telescope's beam. The 
largest pitfall in the virial method is the question of its 
applicability. Mass measurements nearly always neglect 
other contributions that are present in the full virial the- 
orem, such as external pressure, changing moment of in- 
ertia, magnetic fields, the degree of virialization and the 
measurement of a single size for a triaxial system. These 
deviations are frequently parameterized using the virial 
parameter (see above) which is surprisingly constant for 
massive molecular clouds {a ~ 1.5, SRBY, HCS). Thus, 
the virial mass estimate provides a reasonable measure- 
ment of a cloud's dynamical mass. With all these po- 
tential sources of error, the masses of molecular clouds 
are highly uncertain, often to 50%, and this uncertainty 
should be included in the determination of the mass spec- 
trum parameters. 

2.2. Cumulative Mass Spectra 

When a sample contains only a small number of clouds 
[Nciouds < 500), it is still possible to derive the parame- 
ters of a mass spectrum by fitting th e cumulative distri- 
bution of masses. Recent work by iFukui et al.l l|2001f) 
demonstrated the utility of this method for clouds in 
the LMC. The principal difficulty in using this method 
arises in assessing errors to the data in the cumulative 
mass spectrum. Uncertainties appear both in the mass of 
the cloud and in the variance of a random sample being 
drawn from an infinite parent distribution. Practically, 
this results in fitting a truncated power-law function to 
data with errors in both coordinates. The mass coordi- 
nate has an uncertainty from the measurement error and 



the cumulative number has an uncertainty characterized 
by a counting error, equal to \/iV. 

To fit the data, we use the "error-in-variables" method 
for parameter estimation in non-linear functions that 
have uncertai nties in both coo rdinates. The method was 
developed bv iBritt fc Luecke (1973) which, in turn, is 
the full development of a method originally suggested by 
(Peming (1943). An equivalent rnethod was developed 
into an algorithm bv lReillv et al.l l|1993fl which has been 
incorporated into StatLib^. It is this algorithm upon 
which the present work is based although the error-in- 
variables method was also presented to the astronomical 
community with the work of Jeffcrvs ( 1980). The method 
maximizes the likelihood that a set of data (M, N) with 
associated uncertainties {(Jm,o'n) can be drawn from 
a distribution with parameters {Nq, Mo,j}. Since the 
equations of condition cannot be solved algebraically for 
the parameters (as they can in the linear case), the 
minimization is performed iteratively in two interleaved 
phases. First, the true values of the data {i.e. without 
measurement errors) are estimated by maximizing the 
likelihood of being drawn from a distribution with some 
initial guess of parameters. Then, using the estimate of 
the true values of the data, the optimal values of the pa- 
rameters are determined. The process is iterated until 
estimates of both the true data values and the parame- 
ters are determined. 

Instead of performing the fit to the data with the model 
given by equation |21 we use the algebraically equivalent 
expression 

= + 03 (5) 

to improve independent estimates of Mq and A^o which 
are highly covariant in the original formulation. Once 

^ |http : //lib ■ Stat ■ emu ■ edu/ 1 
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the algorithm has converged on a vector of parameters, 
9, we transform the elements of 9 back to the parameters 
of interest. We use a bootstrapping technique to estimate 
uncertainties in the derived parameters, using 100 trials 
to sample the distribution of derived parameters which 
is often non-Gaussian. The quoted values of the uncer- 
tainties in the parameters are the median absolute devia- 
tions of the transformed parameter distribution from the 
bootstrapping trials. Examination of the parameter dis- 
tributions using a large number of bootstrap trials shows 
that the medians adequately characterize the uncertain- 
ties. For some distributions, there are more high-mass 
clouds than expected from the distribution at lower mass 
(i.e. the opposite of a truncation). In this case, the pa- 
rameter 03 converges to zero. When this occurs, we fit a 
power law to the distribution of form 

/ l\f \ 7+1 

N(M' > M) = (6) 

and report only Mq and 7. 

We validated our method by fitting the model to 
random data drawn from power-law distributions with 
known parameters. The trial data have normal deviates 
of known dispersion added to them that simulate the ef- 
fects of measurement error. In these simulations, we find 
that the method both recovers the properties of the dis- 
tribution without bias and produces error estimates from 
bootstrapping that agree well with the scatter in derived 
parameters around the known parameters. This implies 
that we are properly accounting for the error in the sam- 
ple as well as recovering the properties of the underlying 
distribution. These tests demonstrate that the error-in- 
variables fit to the cumulative mass distribution should 
be favored over a fit to the binned mass distribution. 

2.3. Systematic Effects 

In addition to the errors in the mass measurement, 
there are also systematic errors in the generation of mass 
spectra. The two dominant contributions to the system- 
atic errors are choice of the mass measurement (virial 
vs. luminous) and the method used to generate the cloud 
catalog. SRBY report Myr oc M^^j^ in their sample 
which implies that determinations of the index 7 can 
vary by 10% depending on the mass measurement. 

The process used to generate the catalogs is likely 
the dominant systematic in measuring the parameters 
of the mass distribution. In particular, the result- 
ing parameters of mass distributions depend on the 
algorithm which assigns flux into the physically sig- 
nificant substructures for which the masses are deter- 
mined. Such decompositions include (1) human assign- 
ment into clouds (e.g. Wilson fc Scoville .1990i) (2) as- 
signment by grouping neighboring pixels above a cut in 
brightness (SRBY, HCS) and (3) computer algorithms 
such as CLUMPFINDjWiUiams ct al. 1994) or GAUSS- 
CLUMPS ijStutzki fc Glisten 1990). Assigning multiple 
distinct structures into a single cloud artificially drives 
the index of the mass spectrum towards more positive 
values. Such blending is most likely to occur when using 
kinematic data to detangle emission in the inner Milky 
Way. Conversely, overzcalous decomposition of objects 
can erroneously split high mass objects into lower mass 
objects, decreasing the value of the index or creating an 
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Fig. 2. — The mass distribution of the SRBY virial mass mea- 
surements. A truncated power-law fit to the data using the meth- 
ods of this study is shown as a solid line. The data show a signif- 
icant break around TV = 50 and the fit recovers this feature well. 



artificial truncation in the distribution. Predicting the 
quantitative impact of these systematics is beyond the 
scope of this work. 

Ideally, the mass distribution should be derived using 
the same decomposition algorithms and mass determi- 
nations from both observations and simulations to min- 
imize these systematic effects. However, the magnitude 
of these systematic effects can be estimated by analyz- 
ing the same data set with different methods. Using the 
derived parameters from the heterogeneous data sets in 
this study, we find that the index can be robustly deter- 
mined, in spite of these systematic effects. Identifying 
truncations and maximum masses are complicated by 
these systematic effects and require care to accurately 
recover (see below). 

3. LOCAL GROUP MASS SPECTRA 

Using fits to the cumulative mass distribution, we have 
reanalyzed the catalogs of GMCs in the Local Group. 
Our results show significant differences in the mass dis- 
tributions of the CMC populations. For each of the cata- 
logs discussed below, we fit a power-law to the cumulative 
mass distribution, including a truncation if appropriate 
(see above). Unless otherwise stated, when both virial 
and luminous measurements of the mass are reported, 
we use an error equal to half the difference between the 
two mass measurements plus a 10% flux calibration er- 
ror added in quadrature. The results of the new fits to 
the Local Group mass distributions are summarized in 
Table ^ The reported errors are the median absolute 
deviation of the derived parameters for 100 bootstrap- 
ping trials. To illustrate a fit to the da ta, we plot the 
result s of the fit to the virial mass data of iSolomon et alJ 
llT987i SRBY) in Figure H 

3.1. The Inner Milky Way 
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TABLE 1 

Parameters of Mass Spectra for GMCs in Local Group Studies 



Object 


Name 


Type 


Num. 


7 


No 


Mo/(lo5 Mq) 


Inner MW 


SRBY 


VT 


190 


-1.53 ±0.07 


36. ± 12. 


29. ± 5.0 


Inner MW 


SRBY 


CO 


173 


-1.53 ± 0.06 


27. ± 11. 


41. ± 9.5 


Inner MW 


SYSCW 


VT 


107 


-1.58 ±0.15 


14. ± 10. 


26. ± 7.6 


Inner MW 


SYSCW 


CO 


97 


-1.41 ± 0.12 


21. ± 13. 


29. ± 7.2 


Outer MW 


HCS 


VT 


227 


-2.56 ± 0.11 




3.2 ±0.78 


Outer MW 


HCS 


CO 


81 


-2.06 ± 0.15 




6.3 ± 3.1 


Outer MW 


BKP 


VT 


336 


-2.29 ± 0.08 


4.5 ±3.5 


2.9 ± 1.0 


Outer MW 


BKP 


CO 


81 


-2.16 ± 0.17 


2.7 ± 2.9 


2.0 ± 1.0 


M33 


EPRB 


CO 


58 


-2.85 ± 0.36 


2.5 ±2.7 


8.6 ± 3.3 


LMC 


NANTEN 


VT 


44 


-1.71 ± 0.19 


10. ±6.5 


23. ±4.6 


LMC 


NANTEN 


CO 


55 


-1.72 ± 0.12 


6.1 ±3.6 


82. ± 32. 



*The mass distribution shows an excess of clouds at high mass, implying there is no 
truncation in the sample. A pure power law has been fit to the data (Equation ISl. 



There are two major studies of GMCs in the in- 
ner Milky Way. Both SRBY and iScoville et aP l(T987l 
SYSCW) analyzed FCRAO survey data from the first 
quadrant of the Galaxy decomposing the emission into 
clouds using different algorithms. Comparing the results 
of these two studies highlights the systematic effects of 
using different decomposition algorithms. Both studies 
identified clouds as contiguous regions above a fixed an- 
tenna temperature cutoff but chose different thresholds 
and methods for decomposing substructure. We use their 
measurements for virial mass and luminous mass, cor- 
recting for differences in virial definitions and galactic 
scales as summarized in iWilliams fc McKe3 l(l997j) . The 
index of the power-law is unaffected by the choice of con- 
version factor. Typical mass errors are factors of ~ 15%. 
We fit all clouds with masses greater than 1 x 10^ Mq 
in the SRBY study, which approximates their reported 
completeness level. For the SYSCW study, we compared 
the virial and luminous mass measurements after scaling 
the data and we find that the virial mass estimates are a 
factor of 2 higher than the luminous mass measurements 
for the high mass clouds. To place the samples of equal 
footing, we scaled the luminous mass of the clouds by a 
factor of two to bring the mass estimates into agreement. 
We then examined the distributions and established a 
completeness limit of 5 x 10'* Mq based on where the 
distribution departed from a power law on the low mass 
end. Fitting to both the virial and luminous masses for 
both studies finds an index 7 f« —1.5 and a significant 
cutoff with ~ 25 clouds at the cutoff. For all four fits, 
Mq « 3 X 10^ Mq. The relatively small differences be- 
tween the derived parameters despite the different cata- 
log methods suggests that systematic effects are small in 
this case. Since the cataloging methods are conceptually 
similar in the two studies, this result is not surprising. 
The derived value of Nq is slightly higher in the SRBY 
method than SYSCW suggesting there is some influence 
of the catalog method on the cutoff values. 

The mass distribution for the inner Milky Way is shal- 
lower than found for other systems. Two effects may 
bias the results to a shallower index. First, line-of-sight 
blending will make several less massive clouds appear as 
a single, more massive cloud shifting the index to shal- 
lower values. The methods used to generate the SRBY 
and SYSCW catalogs do little to split up blends of emis- 



sion. Second, incorrectly resolving the distance ambi- 
guity will also bias the mass distribution to shallow in- 
dices. If every cloud has the same probability of having 
its distance incorrectly determined, then more low-mass 
clouds at the near distance will be erroneously counted as 
a high-mass clouds at the far distance than the reverse, 
simply because there are more low-mass clouds. This lat- 
ter bias can increase the index of the mass distribution 
by as much as 0.2 for 20% of clouds being assigned to the 
wrong distance. Thus, the index of the mass distribution 
for the Inner Milky Way very likely is steeper than can 
be derived from the current observational data. 

3.2. Outer Milky Way 

The data used for the Outer Milky Way are from 
the F CRA O survey of a section of the second quad- 
rant f Hever et al. 1998) which were subse quently ana- 
lyzed by both HCS and lBrunt etHl ll2(M BKP). HCS 
used a cloud extraction algorithm similar to SRBY, but 
defined cloud properties from the intensity distributions 
slightly differently. In contrast, BKP used a modified 
CLUMPFIND algorithm to identify peaks in the emis- 
sion distribution as the nuclei of distinct clouds. Their 
algorithm extracts roughly ~ 50% more sources than the 
work of HCS. They assign cloud properties to the emis- 
sion distribution in a similar fashion as HCS. Mass errors 
in the HCS study are given as half the difference between 
mass measurements plus a flux error and errors in BKP 
are reported in their study. Since clouds in the outer 
Galaxy with masses smaller than 10'' Mq are not virial- 
ized, we set 10^ Mq as the lower mass limit for the fits to 
these catalogs. Adopting this truncation includes many 
more virial mass measurements than luminous mass mea- 
surements since the virial mass tends to overestimate the 
mass of clouds with A/lum < 10"* Mq. Thus, the lumi- 
nous mass distribution likely represents the underlying 
mass distribution better than the virial mass distribu- 
tion. We also require the kinematic distance to be larger 
than 2 kpc to minimize errors in the distance determi- 
nation. We find that the index of the mass distribution 
is steeper than reported in HCS, which is due to the 
improved fitting methods (7 = —2.1 vs. — 1.8 in HCS). 
Since the luminous mass is likely a better tracer of cloud 
mass, we also perform a fit to the luminous mass data 
alone using a lower limit of 2 x 10'^ Mq and derive an 



6 



Rosolowsky 



index of 7 = -2.05 ± 0.06. 

The catalog of HCS shows more clouds than would be 
expected at high mass given a power-law extrapolation 
from lower masses. Such an excess is not seen in the 
BKP catalog, because of the more aggressive decompo- 
sition algorithm employed in the latter study. Without 
careful analysis of the individual clouds, it is impossi- 
ble to say what represents the true distribution of clouds 
at high mass in the outer Galaxy. Since evidence for a 
cutoff appears in the BKP catalog but not in the HCS 
catalog, comparing these two studies illustrates the sys- 
tematic effects of different catalog methods. There is not 
the strong evidence for a truncation in the outer Galaxy 
that is found for the inner Galaxy data. This is likely be- 
cause there are too few molecular clouds to populate the 
distribution up to the truncation mass. Nonetheless, the 
index of the mass distribution is well-determined and is 
significantly steeper than that found in the inner Milky 
Way. 

3.3. MS3 

M33 is the only spiral galaxy for which a cata- 
log of GMCs exists w ith a known completeness limit 
ijEneargiola et alJl200l EPRB). Since the galaxy is seen 
from an external perspective, blending effects are dra- 
matically reduced compared to Milky Way studies. How- 
ever, there are only 59 clouds above the reported com- 
pleteness limit of 1.5 X IO^Mq, and the clouds have only 
CO masses reported since the individual clouds arc not 
resolved. A follow-up study (Rosolowskv ct al. 2003) 
shows that the virial mass is proportional to the luminous 
mass for GMCs in M33 and that the Mqq /M^t does not 
vary significantly over the galaxy. We estimate the er- 
ror in their measurements as the difference between the 
measured and corrected mass discussed in EPRB plus 
their quoted 25% calibration error in the flux scale of 
the interferometer. The derived value of the mass index 
(7 = —2.9) is very steep. M33 is also the most distant 
galaxy in this reanalysis and observational biases may 
affect the index of the mass distribution. However, the 
potential biases would only make the mass spectrum ap- 
pear shallower than it actually is. In particular, blending 
effects will make several less massive clouds appear as 
a single massive cloud, and underestimates of the com- 
pleteness limit will cause the number of low-mass clouds 
to be underestimated. The influence of either of these ef- 
fects would imply that the mass index is actually steeper 
than what is measured: 7 < —2.9. 

It is likely that the extremely steep slope of the mass 
distribution is actually the tail of a distribution with 
a cutoff mass below the completeness limit of the sur- 
vey. EPRB estimate a characteristic mass between 
3 — 7 X 10'* Mq, which could simply be a cutoff mass 
in a truncated power-law distribution. To illustrate the 
effects of fitting a truncated power-law distribution above 
the cutoff mass, we repeated the analysis of clouds in the 
inner Milky Way restricting the sample to clouds near the 
cutoff mass (M > 2 x lO*" M©). Fitting to the restricted 
sample gives 7 = —2.2 with no evidence of truncation in- 
stead of 7 = —1.5 with a truncation. This supports our 
conjecture that the steep slope of the M33 mass distribu- 
tion can be attributed to fitting a power-law distribution 
above the mass cutoff. 



3.4. Large Magellanic Cloud 

The only other complete survey of GMCs in a galaxy 
was completed using the NANTEN 4-m telescope to ob- 
serve the LMC. lMizuno et all pOOl report the most re- 
cent catalog of GMCs, including 55 resolved GMCs for 
which virial masses can be measured. A subsequent pa- 
per (Fukui et al. 2001j) reports an index of 7 = —1.9 
using CO and virial masses from a currently unavailable 
catalog of more GMCs. All of the resolved clouds have 
masses above the completeness limit of the survey. Us- 
ing the virial masses for the 55 reported cl ouds, we de- 
r ive a mass spectrum index consistent with lFukui et all 
(|2001|) with some evidence of truncation. The index on 
the mass distribution derived from the virial masses is 
likely a lower limit (i.e. 7 > —1.9) because the reported 
virial mass measurements do not account for beam con- 
volution. The error-in- variables fit to the data finds that 
the mass distr ibution is shallow er (7 = — 1.7± 0.2) than 
reported in .Fukui et al.l ( |2001[ ) with some evidence of a 
cutoff. The maximum mass in the LMC is similar to that 
in the inner Milky Way (3 x 10^ Mq), but the value is 
poorly constrained by the limited number of clouds in 
the catalog. 

4. DISCUSSION 

There is a real variation in the mass distribution of 
GMCs across the Local Group with indices ranging from 
7 = —2.9 to —1.5. There are cutoffs at a maximum mass 
of 10^-^ Mq in catalogs from the inner Milky Way and 
the LMC. In general, the differences in the mass distri- 
butions have been unappreciated or trivialized; but they 
are, in fact, significant. In the inner Milky Way, the top- 
heavy mass distribution means that studying the most 
massive clouds encompasses most of the star-formation 
in that part of the Galaxy. In contrast, low mass clouds 
contain a substantial fraction of the molecular mass in 
the outer Milky Way and M33. In systems with bottom- 
heavy mass distributions, the star-forming properties of 
these low mass clouds must be examined to obtain a com- 
plete picture of the star-forming ISM. Using 7 « —1.5 
is appropriate for the inner Milky Way but not for all 
galaxies. 

Since molecular clouds of a give n mass appear to 
be similar across the L ocal Group l|Hever et al.l I2fl0 it 
iRosolowskv et alJl2003D . variation among the mass dis- 
tributions is one of the only distinguishing features 
among molecular cloud populations. Owing to rela- 
tively short molecular cloud lifetimes llBlitz fc ShiJlT98(l 
iLeisawitz et al.llT989l: lYamaguchi et al.ll200H) . molecular 
clouds have little time to increase significantly in mass 
due to cloud collisions and accretion. However, the de- 
struction of molecular clouds by their stellar progeny will 
change their mass through photodissociation and hydro- 
dynamic effects. Observations show that the star for- 
mation rate scales roughly with cloud mass in the Milky 
Way (Mooncv & Solomon 1988). If this is approximately 
correct throughout the Local Group, then differences in 
the mass distribution of molecular clouds are not likely 
to arise from different star formation rates. It seems 
likely that differences observed in the mass distributions 
must be due primarily to the formation mechanism of 
molecular clouds. Since many studies see k to ex plain 
the mass distribution o f mole c ular clouds iKwanlil979.: 
lElmegreen fc Fakarond 119961: iVazauez-Semadeni et alJ 
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ll997HStutzki et alJll998HWada et aIll200(H) . these expla- 
nations must be expanded in scope to encompass the va- 
riety of mass distributions observed in the Local Group. 

It is interesting to note that the environment with the 
steepest index of the mass distribution (M33) is also the 
region which is most gravit ationally stable with respec t 
to gravitational instability l)Martin fc KennicuttI 120011) . 
The behavior might be expected if two mechanisms dom- 
inated the cloud formation process, each producing dif- 
ferent mass distributions and one of the mechanisms was 
regulated by gravitational instability. For example, if 
the molecular clouds that form in spiral arms are more 
massive than those that form in the field, then a steeper 
mass index is expected in M33 where the disk is stable. 
Another possibility is that the galactic environment es- 
tablishes the cutoff mass for the mass distribution. In 
both the inner Milky Way and the LMC where there is 
reasonably clear evidence for a cutoff mass, that mass is 
roughly 3 x 10^ Mq. However, in M33, the characteristic 
mass of molecular clouds must be smaller than the com- 
pleteness limit in the study (1.5 x 10^ Mq) and is likely 
~ 5 X IC* Mq. The outer Milky Way does not appear 
to show a characteristic mass which can be attributed to 
the absence of sufficient molecular material to populate 
the distribution at masses near the cutoff. It remains 
an open question as to what physics would establish the 
characteristic mass in these systems and why the charac- 
teristic mass in M33 would be two orders of magnitude 
less massive than in the Milky Way and the LMC. 

5. CONCLUSIONS 

This study emphasizes the importance of performing a 
uniform analysis to generate mass spectra. Using the 
error-in-variables method of parameter determination, 
we reanalyzed the molecular cloud catalogs for the Local 
Group of galaxies and we report the following conclu- 
sions: 

1) Fits to the cumulative mass distribution using the 
error-in-variables method produce a reliable estimate of 
the parameters of the mass distribution. Bootstrapping 
produces reasonable uncertainties these parameters. The 
adopted method is superior to the standard technique of 
fitting a binned approximation to the differential mass 
spectrum since it is insensitive to bin selection and it 
also accounts for uncertainties in the mass estimate. 

2) There is significant variation in the mass distribu- 
tions of molecular clouds across the Local Group even 
after accounting for systematic effects and biases. Dif- 
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ferences in the method used to catalog the molecular 
emission affect the derived parameters of the mass distri- 
bution. In particular, the presence and magnitude of a 
cutoff in the mass distribution is affected by the decom- 
position algorithm. Unless the cutoff is quite significant 
(as it is in the inner Milky Way), the presence of a trun- 
cation should be regarded with some suspicion. However, 
the index of the mass distribution is far less sensitive to 
the particulars of the mass determination and decom- 
position algorithm, resulting in systematic errors in the 
index 7 of ±0.1. 

3) The mass distribution in the inner Milky Way has 
a measured index of7 = — 1.5±0.1 with good evidence 
for a truncation in the distribution setting a maximum 
mass of 10^'^ Mq. Systematic errors particular to the 
study of the inner Milky Way suggest that the true mass 
distribution may be steeper than this derived value. Us- 
ing 7 « —1.5 is appropriate for the inner Milky Way but 
does not approximate the mass distribution of molecular 
clouds across all galaxies. 

4) The mass distribution of molecular clouds in the 
outer Milky Way is significantly steeper than that found 
in the inner Galaxy. The mass distribution has an index 
of 7 = —2.1 ± 0.2, steeper than previously claimed, and 
shows no evidence of a cutoff at high mass. 

5) The GMCs in M33 show the steepest distribution 
found in this study with no evidence of a cutoff. It is 
possible that the distribution actually has a cutoff below 
the completeness limit of the sample which accounts for 
the derived index. 

6) The LMC has a mass distribution that is steeper 
than that of the inner Milky Way (7lmc = -1.7 ± 0.2) 
but also shows some evidence of a cutoff near 10^-^ Mq 
which was unknown heretofore. An expanded catalog of 
clouds is needed to confirm this result. 
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